Effective field tfieory with a 6'-vacua structure for 2d spin systems 
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We construct a nonlinear sigma (NLct) model description of 2+ Id spin systems, by coupling 
together antiferromagnetic spin chains via interchain exchange terms. Our mapping incorporates 
methods developed recently by ourselves and by Senthil and Fisher, which aim at describing compe- 
tition between antiferromagnetic and valence-bond-solid orders in quantum magnets. The resulting 
2-|-ld 0(4) NLcr model contains a topological 0-term whose vacuum angle 9 varies continuously 
with 5, the bond-alternation strength of the interchain interaction. This implies that the S-vacua 
structure for this NLcr model can be explored by tuning 5 in a suitable 2d spin system, which is 
strongly reminiscent of the situation for 1+1 AF spin chains with bond-alternation. 
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Haldane's semiclassical description of antiferromagetic 
(AF) spin chains in terms of an 0(3) nonlinear sigma 
(NLct) model is a simple and yet powerful format for 
detecting quantum exotica in ld0. In this language 
one exercises special care to keep track of possible "6'- 
term" ^ type Berry phases, which introduces complex- 
number valued weights, and hence nontrivial quantum 
interference into the the path-integral representation of 
the partition function, Z[n] = / X>ne~^^"'^''~^^'>\ Here 
the 0-term action is Se = idQrx, where the winding num- 
ber Qrx = 4^ / drdxri ■ drfi x dxTi probes the global 
topology of the unit Neel vector n's space-time config- 
uration. Hereafter we refer to the coefhcient 9 as the 
vacuum angle. The presence of Sg has profound conse- 
quences, altering nonperturbatively the system's ground 
state properties. In particular, the system undergoes 
a continuous quantum phase transition when the vac- 
uum angle traverses the value 9 = it modulo 27r. Fol- 
lowing the success in discriminating between the spec- 
tral properties of integer and half-integer spin systems 
0, this aproach was applied to bond-alternating AF 
Heisenberg spin chains Q, described by the Hamiltonian 
H = J2i ~ {—^yS)Si ■ Si+i. The vacuum angle was 
found to depend on the bond-alternation strength S as 
9 = 27rS'(l — S). Sweeping S e [—1,1] would therefore 
enable one to probe the entire 6'-vacua structure. Physi- 
cally, this implies the existence of 25" successive quantum 
phase transitions between different valence-bond-solid 
(VBS) states, as later verified by numerical studies^- 
This mapping technique is now incorporated routinely 
to determine the global phase diagram of various Id and 
quasi-ld AF systems. 

The present study is an attempt to construct a NLct 
model exhibiting a nontrivial 0-vacua structure, which 
describes a class of bulk 2d spin systems. This is in part 
motivated by recent interest in exotic phase transitions in 
2dAFs,e.g. between 2d VBS states^- Here we will show 
that certain inter- VBS transitions in 2d can, like in Id, 
indeed be described by a NLct model with a 9-teim, al- 
beit with the target manifold 0(4). (Observe that four is 
the required number of components for constructing the 



winding number Qrxy (defined later), which is the 2-f Id 
analogue of Qrx-) Searches for novel Berry phase terms 
within in the 0(3) NLct model description of 2d AFs were 
carried out in the past 0,110, which showed that singu- 
lar hedgehog events can induce Berry phase factors that 
drive the system towards a VBS groundstate 0,0 ■ An al- 
ternative approach is to depart from the 0(3) description 
where the amplitude of the AF order parameter is kept 
fixed, and to employ instead a composite 5-component 
order parameter explicitly representing the competition 
between the three AF and two (horizontal and vertical) 
VBS components. We found in a previous work starting 
from the 7r-flux mean- field anstaz 9J , that the resulting 
2-l-ld 0(5) NLct model contains a Wess-Zumino (WZ) 
term which reproduces the hedgehog Berry phase factors 
in certain limits. Here we follow this latter route. 

We had observed in ref.p] that "freezing" the fluctu- 
ating VBS amplitude along one of the spatial directions, 
thereby reducing the number of effective components to 
four, will modify our effective 0(5) action into a 2+ld 
0(4) NLct model with a 0-term. This would then sug- 
gest that a novel 0-vacua structure arises in a system 
with a spatially anisotropic AF-VBS competition, a sit- 
uation which might be realized e.g. in a coupled stack of 
Id spin chains with sufficient intrachain AF-VBS com- 
petition. The direct relation to such systems is not to- 
tally apparent though, since we began with a spatially 
isotropic 2d system. For this reason we shall employ and 
extend a method recently suggested by Senthil and Fisher 
[iQ] . also along similar lines, which starts from purely Id 
systems. Upon coupling these system uniformly along a 
second direction, these authors arrive at the 0(4) model 
with 9 = TT. Below we will find that introducing a tunable 
bond- alternation along the transverse direction allows us 
to access arbitrary values of 9, i.e. a 0- vacua very similar 
to the spin chain problem mentioned above. 

Our starting point is the level-1 1+ld SU{2) Wess- 
Zumino- Witten (WZW) model (higher-level cases, i.e. 
higher-S* generalizations will be mentioned later). Here 
we collect features of this model which will shortly be- 
come necessary. Physically, the WZW action simply de- 
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scribes the quantum criticality arising from competing 
AF and VBS orders |3, E • This becomes transpar- 
ent when one represents the S't/(2)-valued field g using 
a unit 4-vector (j)^ (a = 0, 1, 2, 3), viz g = + icj) ■ cr, 
where (p = (f)2, 4>3)- Indeed, it is straightforward to 
show with the aid of nonabehan bosonization techniques, 
that 00 and <p each corresponds in spin language to the 
dimer (VBS) order parameter (—1)* < Si ■ Si+i >, and 
the Neel vector n^x) Expressed as a functional of 

the composite AF-VBS order parameter (j)a, the WZW 
action is an 0(4) NLa model with a WZ term, 



1 

8^ 



drdx [{dr(l)af + {dx<t>a) 



n^ai (1) 



where the WZ functional is 



abed /•! 



dt / dTdx(l)adt(t>bdT(t>cdx(t>d- (2) 



The extra parameter t G [0,1] is a complication common 
to all WZ-type terms (the simplest example being the 
Berry phase for a single spin); here it is used to locally 
frame the three-dimensional "area" on the surface of the 
hypersphere 6*3 swept out by the map (x, r) G S'2 — > 0^ G 
53. We record two properties of r[0„] which we incorpo- 
rate below. First, the change under a slight deformation 
of the configuration (j)a (f'a + can be written as a 
local functional (i.e. one without t): 



5V[cj)a{T,x)\ = 



abed 



dTdx4)a54>bdr4>cdx4'd 



(3) 



Secondly, since r[0c,] is linear in all four components of 
the vector (/)q,, simultaneously flipping any three will in- 
duce a sign change. In particular. 



-0] = -T[4>Q,4>\. 



(4) 



Let us now follow ref.^3 and stack our Id (/>q,(t, x) 
chains along a second {y) spatial direction. We include 
a bond-alternation 5 in the interchain coupling, antici- 
pating its contribution to the 0-term. As detailed below, 
the interchain interaction iSj^ = J drH^ required for our 
purpose of generating a family of 0-vacua has the form, 

Hi_ = - JdxJ2[M^ + {-^r^)N{x,y)-N{x,y + l) 



+ Jx(l - {-iyS)No{x, y)No{x,y + 1)], 



(5) 



with J_L > 0, and Na{x,y) ^ {No{x,y), N{x,y)) = 
iM^-Jj) + a{-l)yio{x,y),{-l)y<f){x,y) + al{x,y)) {a 
is the lattice constant). Here we have reinstated a 
small and rapidly fluctuating component la{x,y) = 
{lo(x,y),l{x,y)), to be eventually integrated out, on top 
of the slowly varying field 4)a{x,y). This intermedi- 
ate step is easily understood when one recalls that for 
the Id AF, the uniform magnetization L (the subdomi- 
nant fiuctuation in that case), playing an analogous role. 



needed to be integrated out to arrive at the final ef- 
fective actional. To preserve the SU{2) symmetry of 
the WZW field g, we need to impose the constraint 
NaNa ^ + N'^ = I, which implies that (pJa = 0. 
The form of H± dictates how the slow and rapid modes of 
the WZW fields on adjacent chains are to align; the slow 
fluctuations consist of a staggered alignment of 4>{x, y) 
fields and a columnar alignment of 4>o{x,y) fields along 
the y-direction poj . This is consistent with our previ- 
ous work, which suggested that an intrinsic competition 
between 2d AF and columnar dimer states exists in the 
vicinity of the 7r-flux state . As for the rapidly varying 
modes, the vector l{x,y) tends to pile up in a ferromag- 
netic fashion, while the component lo corresponds to a 
staggered VBS configuration. This situation, depicted in 
fig.l, would be natural in the presence of a frustrating 
diagonal exchange interaction which is not dominantly 
strong. Later we will return to the issue of relating our 
study to a real spin system. 



FIG. 1: (Figure appended as a GIF file) Slowly and rapidly 
fluctuating degrees of freedom. Crossed lines indicate the 
presence of a frustrating diagonal exchange discussed in the 
text. 

Our goal below is to extract an effective theory for 
(j)a{T,x,y). Readers familiar with Haldane's mapping 
for AF spin chains |^ may find what follows more 
tractable by keeping in mind the set of correspon- 
dences between the Id and 2d cases, displayed in Ta- 
ble 1. In the second entry in the last row, Qrxy = 

TABLE I: Correspondence with Haldane's mapping in Id 



Id ( 0(3) ) 


2d ( 0(4) ) 


n{T,x) + {~lY%L(r,x) . 


<j)c{T,x,y) + a{-l)yic{T,x,y) 


iSijj[n{T, x)] 


r[(pa{T,x,y)] 


iSuj[n] ~ i J dm ■ Sn x drn 


eq.© 


Lj[—n] = —uj[n] 


eq.lgj 


Se = idQrx 


So = iOQrxy 



J drdxdye'^''^'^'^ (j)adr4'bdx<t'cdy<j3d, which is the winding 
number associated with the homotopy 113 (6*3), enters into 
the e-ievTn for the 2+lA 0(4) NLct model. 

We now proceed to the continuum approximation. We 
begin by carrying out the y-summation over the WZ func- 
tionals with the help of eqs.Q and Q, 

^F[0o + a(-l)^Zo,(-l)'^0 + ai] 
y 

= ;^(-l)«F[0„] + '-^e-'''''Y,fd^dxcbalbdrcl>cdxcl)d-{G) 
y y'' 

The alternating series in the last line is converted into an 
integral [Tfij. 

Y^{~l)yT[My)\ dydyT[4>^] - -i^Qrxy (7) 



3 



Next we turn to the interchain term H± which is de- 
composed into two contributions, H±i and H±2i where 



H±i = - J dxJ± N4y)Nc.{y + 1), (8) 



and 



= ' j dxJ^5Y,{~\)y[N^{y)N^(y^\) 



-N{v)-N{v + l)]. 



(9) 



Here we have suppressed the explicit dependence on x 
for brevity. The continuum form of these interactions 
(discarding osciUatory contributions) read 

iJ^i - / dxdy (^a{dy(l)^f + 2J^all\ , (10) 



and 



H±^2 ^ 2J±Sa / dxdylad. 



(11) 



FinaUy we integrate over la, by coUecting terms from 
eqs.®, l(Tn|) and ((TT|l . and completing the square with 
respect to la- Dropping a higher-order derivative term 
and carrying out a suitable rescaling, our effective action 
reads 



Seff[(l)a{T,x,y)] = ~iTr{l-d)Q 



Txy 



-)- d^x 



2g^v 



where the velocity v and the coupling constant g each 
depends on parameters of the original Hamiltonian. The 
main findings here are (1) the dependence of the vacuum 
angle on the bond-alternation parameter, 9 = 7r(l — 5), 
and (2) the factor (1 — 5"^) which enters the coefficient 
for the interchain kinetic energy, both of which coincide 
with known results 3] for the S=l/2 bond alternated spin 
chain. In particular the vacuum angle = n (corre- 
sponding to (5 = 0) lies at the point in parameter space 
where the 'strong' and 'weak' vertical bonds are inter- 
changed. Further analogy with the 0- vacua for the 1-l-ld 
case arises, when one generalizes the foregoing mapping 
to general S, by starting with the level 25* WZW model. 
This gives 9 = 27r5(l — S), an exact reproduction of the 
Id result mentioned earlier. These strong similaritites nat- 
urally lead us to adopt the physical picture established 
within the spin chain context, and associate the 9 = t: 
{mod 27r) points with phase transitions between different 
{vertical) VBS ground states. (In contrast to the Id prob- 
lem, however, the precise nature (including the order of 
transition) of the 2-|-ld 0(4) NLcr model at = tt so 
far does not seem to have been pursued.) In the pres- 
ence of an AF-favoring anisotropy, one can seek support 
for this expectation by viewing the problem in terms of 



the meron gas expansion '5, described (for 9 = n) in 
ref. . Extending this to arbitrary 9 results in an effec- 
tive sine- Gordon type theory of the form 



S„ 



f 9 

[(p\ — j dTdxdy[K{d^ip)^ + Xi cos -cos if 

+A2 cos6'cos2(^ + • • •]. (13) 



In the above action, each harmonic term (oc cos{nip)) 
represent processes associated with single, doubled, and 
tripled merons, etc. When 9 w {S 0), the cos(/7 
term is the most relevant (provided the fugacity ex- 
pansion is valid) and picks out a unique ground state 
value for ip with (depending on 9) either e^^ = 1 or 
e'^'-f = — 1. Resorting to the symmetry analysis of Scnthil 
and Fisher (10|, one finds that this state corresponds to 
a vertical VBS state which simply follows the externally 
imposed bond-alternation pattern. The special role of the 
9 — TT point manifests itself here in the vanishing of the 
cos if term and all other odd harmonics. Here the doubled 
meron term cos2(^ becomes dominant, leading to doubly 
degenerate ground states satisfying e*"^ = ±i[l(lj. Hence, 
tracing the change of 9 through 9 — tt shows that the 
system switches at this point between two vertical VBS 
states, corresponding to e"'' = 1 and e^"^ = —1 (Fig. 2). 



FIG. 2: (Figure appended as a GIF file)The S-vacua for S=l/2 
in the presence of AF-favoring anisotropy. 

We now discuss how the above 0-vacua structure is 
likely to emerge in the context of actual spin systems. 
We know that 2d AFs have a much stronger tendency to- 
wards Neel ordering than in Id, and one would generally 
expect an AF phase to intervene between different VBS 
phases ^4|. As briefly mentioned earlier, a partial clue 
regarding this question lies in the mapping process itself. 
Recalling that our low energy theory resulted from inte- 
grating out the subdominant modes la, it is clear that 
the system consistent with the derivation should have 
such modes as the chief rapid fluctuations. Upon intro- 
ducing a frustrating diagonal exchange term, configura- 
tions such as depicted in the lower half of Fig. 1 would 
start to have significant weight, thus weakening the AF 
dominance. (Similar situations may also arise when pro- 
jecting a three dimensional anisotropic frustrated magnet 
onto an effective model of coupled chains 0|.) If a di- 
rect inter- VBS transition is observed with an appropriate 
amount of frustration -perhaps supplemented with addi- 
tional perturbations, our results suggest that such points 
are described by the strong coupling regime of our ef- 
fective theory at 9 — n. It would also be interesting to 
relate our framework with large-N studies of frustrated 
magnets ll^. Thus, although considerably fragile in na- 
ture compared to the Id case, we believe the 0-vacua of 
the 2-|-ld 0(4) model should have relevance to actual 
frustrated magnets. 
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We expect that various aspects of O-term physics in 
Id will obviously find counterparts in our 2+ Id prob- 
lem. For instance, the O-teim for an open spin chain con- 
tains a boundary contribution which induces a fractional- 
spin edge state|17|. This would correspond to decoupled 
spin-chain-like objects at the upper/lower ends of our 
system. One can also construct a 3d system with a 9- 
vacua structure by repeating the above steps in one di- 
mension higher. The appropriate starting point for this 
purpose would be the 2-|-ld 0(5) model which we derived 
in ref.[12], containing the WZ term 

r'+'[0VBS,0AF] = J\t J dTclxdye'^''''^ 

(t)adt(l)bdr4'cdx4>ddy(j)e, (14) 

where the VBS portion of the five-component composite 
order parameter, (/>vbs now consists of two components. 
Comparing with eq.(|2Il, we see that this is the incarna- 
tion of the WZW term in 2+ Id. Stacking up the 2d sys- 



tems along a third (z-)direction, we obtain in the absence 
of bond-alternation a Berry phase term, inQrxyzi where 

Qrxyz = sfr / dTdxdydze°'^'"^''4>adT(t>bdx4>cdy(t>ddz(j>e- 
This is again a NLcr model with 9 = 7:, this time in 
3+ld and on the target manifold 0(5). It is easy to ver- 
ify that adding on an interplanar bond-alternation would 
shift the value of 9 in this case also, which may describe 
novel physics in 3d magnets. 

In summary we have shown that direct transitions 
among VBS states in anisotropic 2d AF systems can 
in certain cases be described in terms of the 0-vacua 
structure of the 2-|-ld 0(4) NLcr model, in much the 
same way that inter- VBS transitions in spin chains have 
been understood using the 0-term of the 1-t-ld 0(3) NLtr 
model. We suggested that realizations may be found in 
anisotropic frustrated magnets. 
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